AD-A102  872  ROYAL  AIRCRAFT  ESTABLISHMENT  FARNBOROUGH  (ENGLAND)  F/G  9/1 

SIMPLE  METHODS  FOR  CALCULATING  THE  AZIMUTHAL  COVERAGE  OF  HF  DIR— ETC(U> 
FE8  81  G  MAY 

UNCLASSIFIED  RAE-TR-81018  DRIC-TR-81Q18  NL 


LEVELS 


GR7812K 

TR  81018^ 

<? 


ROYAL  AIRCRAFT  ESTABLISHMENT 


Technical  81018 


Febru«np#081 


DTfO 

kELECTE| 

.  AUG  1  4  1981 


Simple  jmethods  for^alculating 

'  THE  ^ZIMUTHAL  jgOVERAGE  OF 

|  jiF  Directional  antennas. 


Cl 


(mi 

xcj 


Procurement  Executive,  Ministry  of  Defence 
Farnborough,  Hants 


Approved  for  public  n 
Distribution  UnSni! 


O 


UDC  621.396.67  :  621.3.029.6  :  526.63  :  514.116.3  :  681.321 


ROYAL  AIRCRAFT  ESTABLISHMENT 


Technical  Report  81018 
Received  for  printing  10  February  1981 


SIMPLE  METHODS  FOR  CALCULATING  THE  AZIMUTHAL 
COVERAGE  OF  HF  DIRECTIONAL  ANTENNAS 

by 

G.  May 


SUMMARY 

\  ■  .  —  —  -  ■ 

'-'Methods  of  calculating  the  azimuthal  coverage  of  HF  directional  antennas  are 
described  and  worked  examples  are  given.  Particular  attention  is  paid  to  spherical 
trigonometry  formulas  employing  auxiliary  quantities  because  of  their  suitability  for  use 
with  simple  hand-held  electronic  calculators.  The  relationships  between  sides  and  angles 
of  oblique  spherical  triangles,  collated  from  several  sources,  are  listed  and  the  value 
of  applying  tests  based  on  these  relationships  prior  to  making  detailed  solutions  of 
spherical  triangle  problems  is  demonstrated.^- 
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1  INTRODUCTION 

With  the  increasing  use  of  directional  antennas  for  HF  communications  on  point-to- 
point  and  ground-air  links  the  determination  of  the  azimuthal  coverage  obtained  at  a 
given  range  and  bearing  from  the  ground  antenna  is  of  considerable  interest  and  import¬ 
ance.  Whilst  methods  of  calculating  the  distance  and  bearing  between  two  locations  on 

1  2 

the  earth's  surface  can  be  found  in  cotanonly  available  mathematical  ’  and  radio 
3  4 

engineers'  handbooks  *  ,  it  is  apparent  that  a  number  of  simple  aids  for  resolving  the 
associated  spherical  triangle  problems  are  not  generally  known.  The  reason  for  this 
appears  to  be  that  these  aids  are  scattered  amongst  a  number  of  texts.  Consequently  an 
attempt  has  been  made  to  gather  these  aids  together  along  with  the  conventional  methods 
and  to  demonstrate  their  application  and  value. 

The  advent  of  the  cheap  electronic  hand-held  scientific  calculator  has  also  meant 
that  several  methods,  which  in  the  past  would  have  been  very  time-consuming  and  tedious 
to  use,  can  now  be  employed  swiftly,  simply  and  accurately.  As  an  example  the  use  of 
the  haversine  formula,  which  previously  would  have  required  access  to  haversine  tables, 
can  now  be  readily  employed  since  the  conversions  between  haversine  and  cosine  functions 
can  be  carried  out  using  a  few  key  strokes  on  the  electronic  calculator.  It  is  of 
interest  to  note  that  one  of  the  greatest  British  mathematicians,  John  Napier"*  (1550-1617) 
Eighth  Laird  of  Merchiston,  near  Edinburgh  not  only  invented  logarithms  but  also  Napier's 
analogies  or  relationships  in  spherical  trigonometry  (see  Appendix  A) .  Moreover  he  was 
the  first  to  attempt  the  construction  of  a  calculating  machine  for  the  performance  of 
multiplication  and  division  by  means  of  metal  plates  in  a  box. 

Basic  formulas  used  in  solving  spherical  triangle  problems  are  stated  in  Appendix  A 
while  in  Appendix  B  useful  relationships  between  the  sides  and  angles  of  oblique  spherical 
triangles  are  given.  In  Appendix  C  the  solutions  of  oblique  spherical  triangle  problems 
using  auxiliary  quantities,  a  method  very  appropriate  for  use  with  electronic  calculators, 
are  derived. 

The  use  of  the  formulas,  detailed  in  Appendices  A  and  C,  in  the  calculation  of  range 
and  bearing  between  two  locations  are  discussed  in  section  2  while  detailed  calculations 
and  examples  are  given  in  Appendix  D,  where  the  usefulness  of  the  relationships  stated  in 
Appendix  B  is  illustrated. 

In  section  3  consideration  is  given  to  the  problems  of  finding  the  azimuthal  cover¬ 
age  obtained  by  the  use  of  directive  HF  antennas  and  this  is  amplified  by  the  worked 
examples  of  Appendix  E. 

The  choice  of  units  for  the  calculations  of  Appendices  D  and  E  is  not  easily  made 
owing  to  the  practical  use  in  military  and  world-wide  airline  operations^  of  both 
Imperial  and  metric  units.  The  convention  has  been  adopted,  therefore,  of  calculating 
ranges  and  distances  in  nautical  miles.  A  table  of  conversion  factors  and  constants  is 
given  in  Appendix  F. 

A  further  factor  in  the  writing  of  this  Report  has  been  the  preparation  of  revised 
material  on  great  circle  bearings  and  distances  and  new  material  on  the  azimuthal 
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coverage  aspects  of  HF  directional  antennas  for  use  in  the  Air  Publication^  dealing  with 
the  general  principles  of  HF  and  MF  antennas  for  ground  stations. 

2  CALCULATION  OF  RANGE  AND  BEARINGS  BETWEEN  TWO  LOCATIONS  WHOSE  CO-ORDINATES  ARE  KNOWN 

Two  quantities  that  commonly  have  to  be  calculated  in  planning  and  deploying  HF 
antennas  are  the  bearing  and  distance  between  two  points.  Moreover  these  quantities  have 
often  to  be  found  before  other  factors,  such  as  the  azimuthal  coverage  given  by  an 
antenna  having  a  known  beamwidth,  can  be  calculated. 

Fig  I  shows  the  two  locations  P  and  Q  both  of  which  will  be  seen  to  lie  on 

meridians  through  the  North  Pole.  Thus  by  joining  P  and  Q  by  the  great  circle  through 

these  two  points,  one  has  only  to  solve  the  spherical  triangle  NPQ  in  order  to  determine 

the  value  of  PQ  and  the  two  angles  NPQ  and  N($P .  The  latter  enable  us  to  find  the 
bearings,  from  true  North,  of  one  point  to  the  other,  care  being  taken  to  remember  that 
by  convention  bearings  are  always  taken  in  a  clockwise  manner  from  North  (see  Appendix  D, 
section  D. 1 ) . 

Several  methods  of  solving  the  problems  of  distance  and  bearing  are  demonstrated  in 
Appendix  D  and  a  choice  can  be  made  depending  upon  whether  one  has  access  to  either  a 
small  pocket  electronic  calculator  or  a  set  of  mathematical  tables. 

It  should  be  noted  that  in  Appendix  D  all  the  calculations  of  great  circle  distances 

and  bearings  are  made  on  the  assumption  that  the  earth  is  truly  spherical.  Whilst,  as  is 

pointed  out,  this  assumption  is  not  strictly  correct  it  is  nevertheless  a  reasonable  one 

for  the  reasons  given.  Very  accurate  calculations  of  the  distance  between  two  widely 

separated  points  on  the  earth's  surface,  ie  the  length  of  the  geodesic  line,  can  be  made 

8 

but  it  is  not  a  simple  matter  and  is  beyond  the  scope  of  this  Report.  Also,  in 
Appendix  D,  section  D.4  the  value  of  using  the  relationships  between  the  sides  and  angles 
of  the  spherical  triangle,  as  listed  in  Appendix  B,  is  demonstrated.  These  relationships 
help  to  reduce  the  uncertainties  as  to  which  quadrant  and  hence  which  sign  of  the  tri¬ 
gonometrical  function  should  be  taken. 

3  CALCULATION  OF  THE  AZIMUTHAL  COVERAGE  OF  HF  DIRECTIVE  ANTENNAS 
3. I  Azimuthal  coverage  parameters 

The  need  to  reduce  the  level  of  unwanted  signals  and  interference  has  led  to  the 

9  jo  1112 

use  of  directional  antennas  such  as  rhombics  ,  log-periodic  and  Beverages  *  ,  both 

for  transmission  and  reception.  From  the  specifications  for  these  directional  antennas  one 
can  obtain  their  azimuthal  beamwidth  in  degrees;  typically  the  half  and  tenth  power  beam- 
widths  are  quoted.  The  problem  then  arises  as  to  the  azimuthal  coverage  that  will  be 
obtained  for  the  stated  performance.  Typical  questions  that  arise  are  as  follows :- 

(a)  Given  the  beamwidth  in  degrees  what  will  be  the  latitude  and  longitude  of 
the  points  at  the  extremities  of  that  beam  at  a  given  range? 

(b)  What  is  the  latitude  of  the  point  where  an  extremity  of  the  beam  crosses  a 
known  meridian? 


V 
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(c)  Wbat  is  the  longitude  of  the  point  where  the  beam  extremity  crosses  a  known 
parallel  of  latitude? 

A  diagram  showing  the  main  features  of  these  azimuthal  coverage  problems  is  shown 
in  Fig  2  and  represents  a  directional  antenna  assumed  to  be  located  at  Exeter,  point  P  . 
The  maximum  response  of  the  main  lobe  is  directed  towards  Bod<f>,  point  Q  ,  in  Northern 
Norway.  The  lobe  is  assumed  to  be  symmetrical  and  to  have  a  total  beamwidth  of 
20  degrees,  ie  ±10°  from  the  maximum  response  along  PQ  . 

Calculations  showing  how  the  questions  posed  in  (a),  (b)  and  (c)  above  may  be 
answered  in  this  particular  case  are  given  in  Appendix  E  and  are  summarized  below. 

3.2  To  find  the  latitude  and  longitude  of  the  points  at  the  extremities  of  the  beam  at 
a  given  range 

In  Fig  2,  J  and  K  are  the  points  where  the  ±10°  beam  extremities  cross  the  line 
at  right  angles  through  the  antenna  centre-line  at  Bodifi,  point  Q  .  It  is  shown  in 
Appendix  E  that  the  range  PQ  from  Exeter  to  Sodip  is  1 155  nautical  miles  and  the 
bearing  NPQ  of  Bod<j>  from  Exeter  is  21.9349°. 

The  latitude  of  point  J  ,  X^  ,  is  found  from  the  co-latitude  90  -  Xj  ,  ie  the 

arc  NJ  .  The  longitude  <j>j  is  obtained  using  the  angle  PNJ  since  this  is  the  angle 

from  the  meridian  through  Exeter,  whose  longitude,  if  ,  is  known. 

P 

The  calculations  of  Appendix  E  show  [equations  (E-3)  and  (E-4)]  that  the 
co-ordinates  of  point  J  are 

\j  =  69°10'N  ,  ifj  =  06°57’E  . 

Similarly  from  (E-7)  and  (E-8)  the  co-ordinates  of  K  are 

XR  =  65°05'N  ,  4>k  =  20°33'E  . 

3.3  To  find  the  latitude  of  the  point  where  the  extremity  of  the  antenna  beam  crosses 
a  known  meridian 

In  Fig  2  let  PJ  be  the  northern  limit  of  the  beam  which  intersects,  at  point  H 
(latitude  X  )  ,  the  meridian  NQ  through  Bodif  .  Here  we  have  to  find  NH  the 

n 

co-latitude  of  H  .  Solving  for  NH  using  the  two  angles  and  included  side  method, 
case  2  of  Appendix  C,  it  is  shown  [see  equation  (E-9)]  that  NH  =  15.8993°  and  thus 
XH  -  90  -  NH  »  74. 1007°N  =  74°06'N  . 

3.4  To  find  the  longitude  of  the  point  where  the  extremity  of  the  antenna  beam  crosses 
a  known  parallel  of  latitude 

The  antenna  beam  PJ  is  assumed  to  cross  the  parallel  of  latitude,  X^  ,  at 
point  D  .  The  intersection  of  the  meridian  NP  and  the  parallel  of  latitude,  X^  ,  is 
at  point  E  (Fig  2). 

The  longitude,  <p^  ,  of  D  can  be  obtained  if  we  find  the  angle  PND  which  is  the 
difference  of  longitudes  of  P  and  D  .  In  Appendix  E  it  is  shown,  equation  (E-14), 
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that 


Can 


PND 


sin(Xp  -  Xp)  tan  DPN 
cos  X„  ~ 


and  the  value  of  PND  is  8.7854°  . 

Hence 

d>D  *  4>p  -  PND  -  4.25°W  -  PND  =  4.5354°E  =  04°32'E  . 

4  CONCLUSIONS 

Various  approaches  to  obtaining  great  circle  distances  and  bearings  have  been 
examined  and  the  advantage  of  the  auxiliary-quantities  method  when  using  an  electronic 
calculator  has  been  shown.  The  methods  of  solving  the  spherical  triangles  that  occur  in 
finding  the  azimuthal  coverage  of  HF  directional  antennas  have  been  discussed  and  the 
value  of  examining  the  relationships  between  the  sides  and  angles  of  the  triangles,  prior 
to  making  detailed  calculations,  has  been  demonstrated. 
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Appendix  A 

FORMULAS  FOR  SOLVING  SPHERICAL  TRIANGLE  PROBLEMS 


Adopting  the  usual  conventions  the  sides  (arcs)  are  denoted  by  a,  b,  c  and  the 
angles  opposite  them  by  a ,  8  and  y  respectively  (see  Fig  1). 

Law  of  sines 

sin  a  _  sin  b  _  sin  c  (A-l) 

sin  a  sin  8  sin  y  '  ' 

Law  of  cosines  for  sides 

cos  a  =  cos  b  cos  c  +  sin  b  sin  c  cos  ct  (A-2) 

cos  b  =  cos  c  cos  a  +  sin  c  sin  a  cos  B  (A-3) 

cos  c  =  cos  a  cos  b  +  sin  a  sin  b  cos  y  .  (A-4) 


Law  of  cosines  for  angles 


cos  a  =  -  cos  8  cos  y  +  sin  8  sin  y  cos  a 


cos  8  «  -  cos  y  cos  a  +  sin  y  sin  a  cos  b 


cos  y 


cos  a  cos  8  +  sin  a  sin  8  cos  c 


(A-5) 

(A-6) 

(A-7) 


Napier's  analogies 


tan 


a  -  8 


; in  | (a  - 
iin  (a  + 


sm  i(a  -  b)  y 

1?)  COt  2 


(A-8) 


a  +  8  cos  J(a  -  b)  y 

tan  — r —  <=  - '  T  it  cot  n 

2  cos  }(a  +  b)  2 


(A-9) 


a  -  b 
tan  — j — 


sin  j(a  -  8)  c 

sin  j(a  +  8)  2 


(A-l 0) 


tan 


a  +  b  cos 


}  (a  ~ 
!  (a  + 


L  tan  £ 

cos  J(o  +  8)  2 


(A-l 1 ) 


or 

o 


other  forms  can  be  obtained  by  cyclic  changing  of  the  letters. 
Gauss’  or  Delambre's  analogies 


sm 


a-8  sin  i(a  -  b)  Y 

- s  - — , - -  COS  — 


sin  c/2  C0S  2 


(A-l 2) 
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Appendix  A 


.  a  +  6 

Sin  - r - 


cos  l (a  -  b)  _  y 


cos  c 


71- 


cos 


a  -  8 


sin  Ha  +  b)  y 


sm  c 


7T 


sin  j 


a  +  8  _  cos  i(a  +  b)  y 

— o —  - -  sm  — 

2  cos  c/2  2 


Again  ocher  forms  can  be  obtained  by  cyclic  changing  of  the  letters. 

1 3 

Haversine  formulas  for  the  sides 


hav  a  =  hav(b  -  c)  +  sin  b  sin  c  hav  a 


hav  b  =  hav(c  -  a)  +  sin  c  sin  a  hav  8 


hav  c  =  hav(a  -  b)  +  sin  a  sin  b  hav  y  , 


where  hav  9  is  the  haversine  of  an  angle  9  and  is  defined  by 


hav  9 


J(1  -  cos  9) 


sm 


9/2 


from  (A- 19) 


and 


Haversine  formulas 

hav  a 

hav  8 

hav  y 

where  s  -  J(a  +  b  +  c)  . 


cos  9  =  1-2  hav  0 

0  =  arc  cos(l  -  2  hav  9)  . 

for  the  angles 

=  sin(s  -  b)  sin(s  -  c)  cosec  b  cosec  c 

=  sin(s  -  c)  sin(s  -  a)  cosec  c  cosec  a 

=  sin(s  -  a)  sin(s  -  b)  cosec  a  cosec  b 


(A- 1  3) 

(A- 14) 

(A-l 5) 


(A- 16) 
(A-l 7) 

(A-l 8) 

(A- 19) 

(A- 20) 

(A-21 ) 

(A-22) 

(A-23) 

(A-24) 
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RELATIONSHIPS  BETWEEN  THE  SIDES  AND  ANGLES  OF  OBLIQUE  SPHERICAL  TRIANGLES 

In  solving  problems  associated  with  oblique  spherical  triangles  the  following 
relationships ^  between  the  sides  and  angles  are  found  to  be  of  considerable  value. 
It  should  be  noted  that  they  are  restricted  to  those  triangles  in  which  each  side  and 
angle  is  less  than  180°. 

Bl  The  sum  of  any  two  sides  is  greater  than  the  third. 

B2  The  sum  of  the  three  sides  is  less  than  360°. 

B3  If  two  sides  are  equal,  the  angles  opposite  are  equal  and  conversely. 

B4  The  sum  of  the  three  angles  is  greater  than  180°  and  less  than  540°. 

B5  Half  the  sum  of  two  sides  and  half  the  sum  of  the  two  angles  opposite  them  belong 
to  the  same  quadrant. 

B6  If  a  side  (or  an  angle)  differs  from  90°  by 
side  (or  angle)  then  the  first  side  and  its 
quadrant . 


a  larger  number  of  degrees  than  another 
opposite  angle  belong  to  the  same 


L 


i 
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Appendix  C 

SOLUTION  OF  OBLIQUE  SPHERICAL  TRIANGLE  PROBLEMS  USING  AUXILIARY  QUANTITIES 

In  many  problems  it  is  required  to  find  only  one  or  two  of  the  unknown  parts  of 

a  spherical  triangle.  It  is  evident  from  Appendix  A  that  any  part  can  be  found  using 

any  three  other  parts  through  the  use  of  the  general  formulas  given  there.  By  the 
.  1 4 

introduction  of  auxiliary  quantities  the  basic  sine  and  cosine  formulas  have  been 
adapted  to  logarithmic  computation  and  are  therefore  in  a  form  that  can  be  easily 
processed  on  an  electronic  calculator.  As  these  methods  may  not  be  generally  known 
a  full  statement  will  be  given  of  the  derivation  for  one  case,  two  sides  and  the 
included  angle.  Formulas  for  the  remaining  cases  which  can  be  derived  in  a  similar 
fashion  are  also  given.  Again  the  usual  conventions  are  employed  (see  Appendix  A  and 
Fig  I). 

Case  I .  Two  sides  and  the  included  angle 

Let  a,  b  and  y  be  the  given  parts. 


To  find  c  from 

(A-4)  we  have  that  the 

relation  between  a. 

b,  y  and  c  is 

cos  c  =  cos  a  cos 

b  +  sin  a  sin  b  cos 

Y  .  (C-l) 

Let 

m  sin  M  = 

sin  b  cos  y 

(C-2) 

and 

m  cos  M  = 

cos  b 

(C-3) 

from 

(C-l),  (C— 2)  and 

(C-3)  by  eliminating  b 

cos  c  = 

m  cos  (a  -  M)  . 

(C-4) 

From 

(C-2)  and  (C-3) 

tan  M  = 

tan  b  cos  y 

(C-5) 

and  from  (C-3)  and  (C-4) 


cos  c 


cos  b  cos(a  -  M) 
cos  M 


(C-6) 


To  find  B  Similarly  using  the  relationship  between  a,  b,  y  and  8  it  can 
14 

be  shown  that 


cot  B 


cot  b  sin  a  -  cos  a  cos  y 
sin  y 


and  again  using  (C-5)  it  can  also  be  shown  that 


tan  B 


tan  y  sin  M 
sin (a  -  M) 


(  C—  7 ) 
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To  find  a  Interchanging  a  and  b  and  consequently  a  and  ti  in  ( C— 5 )  and 
(C-7)  and  calling  the  auxiliary  angle  N  we  have 

tan  N  =  tan  a  cos  y  (ti¬ 


tan  a  = 


tan  y  sin  N 
sin(b  -  N) 


cos  c  also  can  be  found  using  N  as  an  auxiliary  angle 


cos  c  = 


cos  a  cos(b  -  N) 


(C-10) 


Case  2.  Given  two  angles  and  the  included  side 


Let  a,  8  and  c  be  the  given  parts. 


To  find  y  and  a  we  have 


tan  a  cos  c 


(C-l 1) 


cos  a  sin(6  -  M) 
sin  M 


(C-l 2) 


tan  c  cos  M 
cos(8  -  M) 


(C-l 3) 


To  find  b  we  have 


tan  N  = 


tan  6  cos  c 


(C-U) 


tan  b  = 


tan  c  cos  N 
cos(a  -  N) 


(C- 15) 


Case  3.  Given  two  sides  and  an  angle  opposite  one  of  them 


Let  a,  b  and  a  be  the  given  parts. 


To  find  c  using 


tan  M  =  tan  b  cos  a 


(C-l 6) 


cos(c  -  M) 


cos  a  cos  M 
cos  b 


(C-l 7) 


To  find  y  and  8  using 


=  tan  a  cos  b 


(C-l 8) 
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and 


sin(y  +  N) 


sin  N  tan  b 
tan  a 


sin  6 


sin  b  sin  a 
sin  a 


(C- 19) 


(020) 


Case  4.  Given  two  angles  and  a  side  opposite  one  of  them 


Let  a,  8  and  a  be  the  given  parts. 


To  find  y 


then 


To  find  c 


then 


tan  M 


sin(y  -  M) 


I _ 

tan  8  cos  a 


cos  a  sin  M 
cos  8 


tan  N 

sin(C  -  N) 


tan  a  cos  8 

tan  8  sin  N 
tan  a 


To  find  b 


sin  b 


sin  8  sin  a 
sm  a 


(021 ) 


(022) 


(023) 

(024) 


(025) 


RIO 


018 
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Appendix  l) 

DETAILED  CALCULATION'S  OF  GREAT  CIRCLE  DISTANCES  AND  BEARINGS  BETWEEN  TWO  POINTS  WHOSE 

CO-ORDINATES  ARE  GIVEN 

D. I  Ceneral  principles 

It  is  desired  to  determine  the  great  circle  distance  and  bearing  between  the  two 
locations  P  and  Q  depicted  in  Fig  1.  It  will  be  seen  that  the  great  circle  arcs  NP 
and  NQ  are  part  of  the  meridians  through  P  and  Q  .  Let  the  latitudes  of  P  and  Q 
be  denoted  by  Xp,  A  and  the  longitudes  by  bp  and  bq  .  In  the  triangle  NPQ, 

NP  =  b  =  90  ±Ap  and  NQ  =  a  =  90  ±Xg  depending  whether  P  and  Q  are  south  or  north 

of  the  equator.  The  angle  PNQ  designated  y  ,  is  derived  from  the  difference  in  longi¬ 
tudes  of  P  and  Q  .  However,  one  must  always  take  the  minor  arc,  thus  y  “  bp  “  9q  or 
360  -  (<)>p  -  bg)  depending  upon  whether  (bp  "  bq)  is  less  or  greater  than  180°. 

To  find  the  great  circle  distance  between  P  and  Q  we  have  to  calculate  the 

value  of  the  arc  PQ  in  degrees,  ie  c  ,  and  then  multiply  by  60  to  obtain  the  answer 
in  nautical  miles.  This  can  also  be  stated  in  kilometres  using  the  internationally  agreed 
conversion  factor'"’  of  1  NM  =  1.852  km.  The  assumption  is  made  in  using  these  conver¬ 
sions,  that  the  earth  is  truly  spherical.  Whilst  it  has  recently  been  determined  that  the 
earth  is  slightly  pear-shaped' ^  the  variations  are  relatively  small  (see  Appendix  E)  and 
are  usually  neglected. 

The  value  of  PQ  can  be  determined  since  we  know  the  values  of  NP,  NQ  (b  and  a 
respectively)  and  the  included  angle  PNQ  ( ie  y) .  Similarly  angle  NPQ  designated  a  , 
or  NQP  designated  8  may  be  determined.  By  convention  true  relative  bearings  are 
always  taken  in  a  clockwise  direction  from  North.  Hence  though  6  is  the  true  bearing 
of  P  relative  to  Q  ,  the  true  bearing  of  Q  relative  to  P  is  360  -  o  . 

D.2  Use  of  cosine  and  sine  laws 

From  equation  (A-4)  we  have 

cos  c  =  cos  a  cos  b  +  sin  a  sin  b  cos  y 


but  a  =  90  ±Aq  ,  b  =  90  ±Xp  ,  y  =  bp  -  bq  or  360  -  (bp  -  bg)  •  When  Xp  and  Xq 
are  both  north  or  both  south  of  the  equator  then 


When  Xp 


cos  c 


is  north  and  X^ 


=  sin  Xp  sin  Xq  +  cos  Xp  cos  Xq  cos  y  . 
is  south  of  the  equator  or  vice  versa 
=  -  (sin  Xp  sin  Xq)  +  cos  Xp  cos  Xg  cos  y  . 


As  an  example  find  the  distance  from  Seattle 

Seattle  (point  P) 
Sydney  (point  Q) 


,  USA  to  Sydney,  Australia 


Lati tude 
47°27'N 

33°54'S 


Longitude 
1  22°21 'W 

1 5 1 0 ) 2'E 


(D-l) 


(0-2) 
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The  difference  in  longitude  is  273°33'  but  this  is  the  major  arc  whereas  as  noted  in 
section  D.l,  y  is  taken  as  the  value  of  the  minor  arc  so 

y  -  360°  -  273°33'  -  86°27'  . 

Substituting  for  X^  ,  X^  and  y  in  equation  (D-l)  and  evaluating  we  find 

cos  c  -  -  0.3761  and  c  -  112.0941° 

thus  the  distance  PQ  is  6726  NM. 

D.3  Use  of  haversine  formula 

A  great  advantage  of  the  haversine  method  is  that  there  is  a  single  positive 
angle  less  than  180°  having  a  given  haversine  value,  thus  eliminating  the  uncertainty 
as  to  which  quadrant  an  angle  or  side  belongs. 

Using  equation  (A-18)  and  decimalising  the  values  for  Seattle  and  Sydney  quoted 
in  section  D.2 


a  -  123.9°  ,  b  -  42.55°  ,  y  -  86.45° 


hav  C  »  hav(a  -  b)  +  sin  a  sin  b  hav  y 

-  hav  81.35°  +  sin  123.9°  sin  42.55°  hav  86.45°.  (D-3) 

Using  (A- 19)  to  find  the  haversine  values 

hav  81.35°  =  0.424  801 
hav  86.45°  =  0.469  040. 

Thus  (D-3)  becomes 

hav  c  =  0.424  801  +  0.263  264 
=  0.688  065  . 

Using  equation  (A-20) 

cos  c  ■  -  0. 376  1 30 
c  -  112.094  15° 

and  the  distance  is  6726  NM. 

D.4  Use  of  auxiliary  quantities  to  find  the  distance  and  bearing  of  one  location  from 
the  other 

In  many  applications  one  has  to  find  only  the  bearing  from  one's  location  to  some 
other  location,  and  it  is  here  that  the  method  described  in  Appendix  C  where  an  auxiliary 
angle  is  used,  proves  to  be  especially  useful. 

Frequently  one  has  to  find  both  the  distance  and  the  bearing  between  two  points  but 
still  only  one  auxiliary  angle  need  be  calculated.  Before  proceeding  with  the  calculation 


it 


Appendix  D 


15 


it  is  advisable  to  ascertain  which  quadrant  the  required  angle  will  be  in  thus  avoiding 
the  uncertainties  that  can  arise.  The  relationships  between  sides  and  angles  of  an 
oblique  spherical  triangle  are  stated  in  Appendix  B  and  an  application  of  these  will  now 
be  shown  using  the  Seattle-Sydney  co-ordinates  again. 

We  know  from  (D-l)  that  we  have  to  find  the  third  side  c  and  either  angle  a  or 
angle  8  (Fig  I)  if  we  wish  to  find  the  distance  and  the  bearing  between  the  two  points 
P  and  Q  .  We  can,  therefore,  use  equation  (C-8)  to  find  the  auxiliary  angle  N  and 
equations  (09)  and  (010)  to  find  both  the  bearing  and  distance,  respectively,  of 
Sydney  from  Seattle  knowing  that 

Q  =  123.90°  ,  b  =  42.55°  ,  y  =  86.45°  . 

Applying  the  tests  in  Appendix  B  we  can  deduce  the  following:-  from  B6  we  see  that 

90  -  b  >  a  -  90  but  b  is  in  the  first  quadrant  hence  8  also  is  in  the  first  quadrant. 

It  then  follows  that  j(f?  +  y)  <  90°  so  from  B5  $(b  +  c)  <  90  and  thus  c  <  137.45°. 

From  B1  (b  +  c)  >  a,  ie  c  >  a  -  b  thus  c  >  81.35°  consequently  81.35°<c<  137.45°. 
We  can  now  see  that  J(a  +  c)  must  be  in  the  second  quadrant  and  hence  again  from  B5 
J(<*  +  y)  is  also  in  the  second  quadrant  thus  180°  <  (a  +  y)  <  360°. 

Now 

Y  =  86.45°  so  93.55°  <  a  <  273.55°  . 

Knowing  also  that  b  and  8  are  in  the  first  quadrant  and  applying  the  sine  rule  we 
see  that  in  order  that  sin  a  may  remain  positive,  a  must  be  in  the  second  quadrant. 
Since  a  and  a  are  in  the  same  quadrant,  the  second,  then  a  -  90  >  c  -  90 
thus  a  >  c  ,  ie  c  <  123.9°. 

Summarizing  we  have 

93.55°  <  a  <  180° 

81.35°  <  c  <  123.9°  . 

Now  from  (C-8) 

tan  N  *  tan  a  cos  y 

*  tan  123.90°  cos  86.45° 

and 

N  -  -  5.2647° 


substituting  in  (C-10) 


and 


cos  c 


cos  a  cos(b  -  N) 
cos  N 

cos  123.90°  008(42.55°  +  5.2647°) 
cos (—  5.2647*  ) 


c  -  112.0941°  . 
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From  (09) 


tan  a 


tan  y  sin  N 
sin(b  -  N) 


and  substituting  we  have 


tan  86.45°  sin(-  5.2647°) 
sin(42.55a+  5.2647°) 


a  =  116.6103°  ■»  1 1 6° 37'  . 


Thus  as  noted  in  (D-l)  the  true  bearing  of  Sydney  from  Seattle,  by  convention  taken  as 
clockwise  from  North,  is 


8(0 
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DETAILED  CALCULATIONS  OF  AZIMUTHAL  COVERAGE  OF  AN  ANTENNA  AT  EXETER  DIRECTED  TO  BODll> 

E . 1  Geometry  of  the  situation 

The  geometry  of  the  situation  is  depicted  in  Fig  2  where  Exeter  is  at  point  P  and 
Bod*  is  at  point  Q  .  The  latitudes  and  longitudes  of  Exeter  and  Bod*  are  as  follows:- 

Latitude  Longitude 

Exeter  50.40°N  4.25°W 

Bod*  67.27°N  14.33°E  . 

Using  one  of  the  methods  described  in  Appendices  C  or  D  the  distance  along  the  great 
circle  PQ  is  found  to  be  19.2431°,  ie  1155  NM  and  the  bearing  of  Bod*  from  Exeter, 

NPQ  ,  21.9349°  from  true  North. 

The  half-power  beamwidth  is  assumed  to  be  20°  thus  JPQ  *  KPQ  =  10° 
where  points  J  and  K  are  on  the  line  JK  at  right  angles  to  the  great  circle  PQ 
and  passing  through  Bod*,  point  Q  . 

Thus 

JOP  *=  KQP  =  90°  . 

Since  NPQ  =  21.9349°,  NPJ  -  11.9349°  and  NPK.  =  31.9349°. 

Also  NP  -  co-latitude  of  P  =  90°  -  50.40°  =  39.60°. 

E.2  To  find  the  latitude  and  longitude  of  points  at  the  extremities  of  the  beam  at 

a  given  range 

For  this  example  we  will  take  the  points  J  and  K  (see  Fig  2) . 

(a)  Point  J 

In  the  spherical  triangle  JPQ  we  know  the  values  of  JPQ,  JQP  and  PQ  so 
we  can  now  use  equations  (C— 14)  and  (C-15).  Noting  that  auxiliary  angle  N  •  0°  since 
tan  JQP  «  tan  90°  =  ”  ,  we  find  tan  PJ  =  tan  PQ/cos  JPQ  . 

This  expression  could  also  have  been  found  from  the  tables  of  formulas  for  right 
spherical  triangles1 ,2’^’ 1 

Thus 

tan  PJ  ■  tan  19.2431°/cos  10° 
so 

PJ  -  19.5177°  ,  ia  1171  NM  .  (E-l) 

Also 

NPJ  -  11.9349°  .  (E-2) 

In  triangle  NPJ  we  now  know  the  values  of  NP,  PJ  and  N$J  ,  thus  using 
equation  (C— 5)  we  find  the  auxiliary  angle  M  to  be  19.1267°  and  we  solve,  using 
equations  (C-6)  and  (C-7)  to  find  NJ  and  PNJ  . 


L 
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From  equation  (D-6)  NJ  =  20.8369°. 

From  equation  (D-7)  PNJ  =  11.2001°. 

Now  NJ  is  the  co-latitude  of  J  and  thus  the  latitude  of  J  is 

=  90°  -  NJ  =  69. !63l°N  =  69°10'N  .  (E-3) 

The  longitude  of  J,  4>j  ,  is  11.2001°  east  of  the  meridian  through  Exeter  which  is  at 
04. 25°W,  ie 

4>j  =  4.25°W  -  11.2001°  =  6 .9501 °E  =  06°57’E  .  (E-4) 

(b)  Point  K 

PK  may  be  calculated  in  the  same  manner  as  for  PJ  .  However,  if  PJ  has  been 
calculated  there  is  strictly  no  need  to  calculate  PK  as  well,  since  by  symmetry  they 
are  equal ,  thus 

PK  =  19.5177°  ,  ie  1171  NM  .  (E-5) 

In  triangle  NPK  we  know,  therefore,  the  values  of  NP,  PK  and  NpK  and  using 
equation  (C-5)  we  find  the  auxiliary  angle,  M  , 

M  =  19.1267°  .  (E-6) 

From  equation  (C-6) 

NK  =  24.9103° 

XR  =  90  -  NK  =  65 .0897°N  =  65°05'N  .  (E-7) 

From  equation  (C-7) 

PNK  =  24.8075°  . 

Therefore  the  longitude  of  K  ,  ,  is  24.8075°  east  of  the  meridian  through  Exeter, 

ie 

<J>K  -  04.25°W  -  24.8075°  =  20.5575°E  =  20°33'E  .  (E-8) 

E.3  To  find  the  latitude  of  the  point  where  the  extremity  of  the  antenna  beam  crosses 
a  known  meridian 

Here  we  will  consider  the  case  where  the  northern  extremity  of  the  beam  crosses 
at  H  ,  the  meridian  through  Bod$.  To  obtain  the  latitude,  X  ,  of  point  H  we  need  to 

H 

find  the  co-latitude  NH  . 

In  the  spherical  triangle  HPN  we  know  the  values  of  the  two  angles  HPN  and 
HNP  and  the  included  side  NP  since  H$N  «  NPJ  «  1 1 .9349°  see  equation  (E— 2)  and 
H#P  ■  PNQ  *  $p  -  $p  ■  18.58°  also  NP  the  co-latitude  of  P  »  39.60°.  Thus  we  can 
calculate  the  value  of  NH  . 
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Using  the  auxiliary  angle  method  of  Appendix  C,  equation  (C— II) 


tan  M 


tan  HPN  cos  NP  tan  1 1 .9349°  cos  39° 


and 


M  =  80.7498  . 


From  (C-13) 


and 


tan  NP  cos  M 

tan  NH  =  - g - 

cos(HNP  -  M) 

NH  =  15.8993°  . 


.*.  -  90°  -  NH  =  74. 1007°N  =  74°06'N  .  (E-9) 

If  we  calculate  the  values  of  NHP  and  HP  using  equations  (C-12),  (C— 14),  (C-15)  but 
first  apply  test  B4  from  Appendix  B  which  shows  that  NHP  >  149. 485°,  we  find 

NHP  =  151.2371°  and 

Applying  the  sine  rule  as  a  check  we  obtain 

sin  15.8993°  sin  39.60° 

sin  11.93490  ~  sin  151.2371° 


HP 


24 .9663 


sin  24.9663 
sin  18.58° 


-  1.324693 


thus  confirming  the  validity  of  the  calculations. 

E .4  To  find  the  longitude  of  the  point  where  the  extremity  of  the  antenna  beam 
crosses  a  known  parallel  of  latitude 

The  northern  extremity  of  the  antenna  beam  PJ  will  again  be  used  and  its  point 
of  intersection,  0  ,  will  be  assumed  to  be  on  the  67th  parallel  of  latitude.  Denote  the 
point  where  the  67th  parallel  crosses  the  meridian,  4>p  ,  through  Exeter,  by  E  .  If 
we  then  find  END  =  PND  and  noting  that  Exeter  is  west  of  the  Greenwich  meridian,  it 
will  be  seen  that  the  longitude  of  point  D  is  given  by 

♦D  -  *p  -  PND  .  (E-10: 

From  the  sine  formula  (A-l)  when  y  *  90°  we  obtain 


sin  a 


sin  a 
sin  c 


and  sin  8 


sin  b 
sin  c 


also  from  the  cosine  formula  (A-4)  when  y 


90° 


COS  c 


cos  a  cos  b 
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combining  the  above  expressions  we  can  show  that 


sin  a  =  tan  b  cot  8  . 


(E-l 1) 


Then  in  triangle  DEN  where  DEN  *=  90  and  EN  =  90  -  X  ,  applying  (E-ll) 


sin(90  -  Xp)  =  tan  ED  cot  DNE 


and  rearranging. 


cos  Xp  tan  DNE  =  tan  ED  . 


(E-l  2) 


Similarly  in  triangle  DEP  , 

Again  using  (E-ll), 


DEP  =  90  and  EP  =  Xp  -  Xp 


sin(Xp  -  Xp)  tan  DPE  =  tan  ED  , 


(E-l  3) 


equating  (E-l 2)  and  (E-l 3)  and  substituting  DPN  for  DPE  since  they  are  identical, 
we  have 

sin(X  -  X  )  tan  DPN 


tan  PND  =  tan  DfiE  = 


D  P' 


cos  X„ 


(E-l 4) 


Now 


67°  ,  Xp  =  50.4°  and  DPN  =  11.935°  . 


tan  PND 


sin  16.6°  tan  11 .935° 


cos  67° 


and 


PND  =  8.7854 


substituting  (E-15)  in  (E-10) 


D  ■  <t>p  -  PND  =  04 ,25°W  -  8.7854°  =  4.5354°E 


04  32'E  . 


(E-15) 


(E-16) 


018 


21 


17 

F.2  Constants 

Earth's  equatorial  radius 
polar  radius 

mean  radius 


6378.160  km  *  3443.845  MM 
6356.775  km  -  3432.384  NM 
6371.020  km  «  3440.076  NM. 
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